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ABSTRACT 

High-redshift quasars (z > 6) drive ionization fronts into the intergalactic medium 
(IGM) , with the thickness of the front generally increasing with the hardness of ion- 
izing spectrum. If the thickness of the front can be measured, it can provide a novel 
constraint on the ionizing spectral energy distribution (SED), uniquely available for 
sources prior to the end of reionization. Here we follow the propagation of an I-front 
into a uniform IGM, and compute its thickness for a range of possible quasar spectra 
and ages, and IGM neutral hydrogen densities and clumping factors. We also explore 
the effects of uniform and non-uniform ionizing backgrounds. We find that even for 
hard spectra, the fronts are initially thin, with a thickness much smaller than the mean 
free path of ionizing photons. The front gradually thickens as it approaches equilibrium 
in 10 8 — 10 9 years, and the thickness of its outer part can eventually significantly exceed 
simple estimates based on the mean free path. With a sufficiently high intrinsic hydro- 
gen column density obscuring the source (log(A^H/cm -2 ) > 19.2) or a sufficiently hard 
power-law spectrum combined with some obscuration (e.g. dlogi^/dlog^^ — 1.2 at 
log(7Vn/cm~ 2 ) > 18.0), the outer thickness of the front exceeds ~ 1 physical Mpc and 
may be measurable from the three-dimensional morphology of its redshifted 21cm 
signal. We find that the highly ionized inner part of the front, which may be probed 
by Lyman line absorption spectra, remains sharp for bright quasars unless a large 
obscuring column (log(-/Vu/cm _2 ) > 19.2) removes most of their ionizing photons up 

to k. 40 eV. For obscured sources with log(./VH/cm -2 ) > 19.8, embedded in a signifi- 
cantly neutral IGM, the black Lyman a trough (where the neutral fraction is ~ 10~ 3 ) 
underestimates the size of the H n region by a factor of > 4. 

Key words: quasars: general - cosmology: theory - observation - diffuse radiation - 
radio lines: general - ultraviolet: general - quasars: absorption lines 



1 INTRODUCTION 

The spectra of z > 6 quasars (e.g. iFan et all 120021 : 
IWhite et al.ll2003l ). showing complete Gunn-Peterson (GP) 
absorption, have been widely analyzed as a possible probe 
of the ionization state of the ambient IGM. The sizes of 



the neutral fraction (ICen & Haiman 2000; Madau & Rees 


2000) and suggest a value above a few percent ( 


Fan et al. 


2002: Wvithe & Loebll2004 iMesinser & Haimanl 


2004), al- 



though this inference is subject to some bias and signif- 
icant scatter, depending on how the intrinsic H II region 
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size i s deduced from the observe d spectrum (|Maselli et al.l 
l2007t iBolton fc Haehneltl |2007al lbh. Direct evidence for a 
GP dampin g wing in the QSO spectra ( Mcs inger fc Haimanl 
2004, 2007) also suggests that the neutral fraction exceeds 
several percent at least along two lines of sight. The evolu- 
tion of Lyman a emitting galaxi es may also sug gest a sig- 
nificant global neutral fraction (lOta et alj 12007 , although 
the presently observed evolution may still be attributed to 
the evol ution in the halo ab undance and in the mean IGM 
density; iDiikstra et""al1l2007l ). 

If the IGM is significantly neutral, then the shape of 
the I-front will be well-defined, i.e. the ionizing flux, over a 
region corresponding to the mean-free-path of ionizing pho- 
tons, will be dominated by the quasar itself, rather than by 
the background. In this case, the shape of the I-front will 
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contain information on the spectrum of the ionizing source. 
In particular, the thickness of the front is generally expected 
to increase for harder source spectra, whose typical ioniz- 
ing photon has a longer m ean free path (hereafter m.f.p.). 
iMesinger fc Haimanl (|2004l ) noted that the onset of the Ly- 
man a and j3 GP troughs in the spectra of z > 6 quasars 
occur at redshifts close to one another, corresponding to a 
physical separation of Rp — R a < 1 Mpc along the line of 
sight, possibly placing an upper limit on the hardness of 
the spectrum (with mean photon energy (E) < 230 eV, for 
which the m.f.p. in a neutral IGM at z ~ 6 is ~ 1 Mpc). Al- 
though the ratio of oscillator strengths of the Lyman a and 
(3 lines is a factor of 5.2, once density inhomogeneities and 
foreground Lyman a absorption is taken into account, the 
Lyman f3 line effectively probes a factor o f 2 — 3 lower neutral 
fract ion than Lyman a (|Fan et al.ll2002l ; ISongaila fc Cowiel 
2002). The comparison of Rp and R a gives the path-length 
over which the neutral fraction rises by a factor of 2 — 3; this 
could naively be taken as a proxy for the thickness of the 
I-front. 

This simple interpretation, however, is complicated by 
possible absorption by the GP dam ping wing from the IGM 
[|Mesinger fc Haimanl |2004 120071 ), and also by the scat- 
ter in the Rp/R a ratio expected to ari se between different 
lines o f sight due to density fluctuations l|Bolton fc HaehnehJ 
2007a). Nevertheless, if a sufficient number of distant 
quasars are detected in the future, the ratio Rp/R a , and 
other features in their absorption spectra, could provide a di- 
agnostic for a finite I-front thickness and therefore the hard- 
ness of the sources, at least statistically. Note that this test is 
not available at lower redshift, where the IGM is highly ion- 
ized, the m.f.p. of ionizing photons is already > 100 Mpc, 
and photoionization equilibrium is established on a time- 
scale of (r/rmAmfp) -1 ~ 10 4 years, much less than the ex- 
pected quasar lifetime. 

Other methods have been proposed to probe th e fi- 
nite thickness of quasar I-fronts. IZaroubi fc Silk (2005) ex- 
plored measuring the thickness of I-fronts through their 
three-dimensional redshifted 21cm signatures, with a low- 
frequency radio telescope array such as LOFAR (for a com- 
prehensive revie w of current and plann ed high-redshift 21- 
cm projects, see iFurlanetto et al.l|200 rj). The soft spectrum 
of stellar radiation would produce a much sharper ioniza- 
tion front than the very hard spectrum predicted to emerge 
from miniquasars; when the sources themselves are not 
detected, the 21cm maps could still di scriminate between 
stellar- and quasar-driven reionization Finally, 
ICantalupo et all (|2007T l recently suggested that collisionally 
excited Lya emission may be detectable from the shell of 
material at the I-front, over an extended solid angle be- 
hind bright quasars; the spectral shape and thickness of this 
emission is another possible measure of the thickness of the 
I-front. 

At present, little is known empirically about the ioniz- 
ing spectrum of quasars at z > 6. At lower redshift, quasars 
are found to have relatively soft power-law spectra in the ex- 
treme ult raviolet, with an a verage spectral index of around 
s = 1.8 (|Telfer et al.ll2002l ). However, many quasars show 
a strong soft X-ray excess at a few lOOeV, which can not 
be explained by thermal emission from a standard accre- 
tion disk. The excess X-rays could originate from Compton 
scattering in a hot corona surrounding the accretion disk 



around the AGN (|Porquet et al.ll2004l ) . Quasars at high red- 
shift are also expected to be preferentially obscured due to 
the increased gas density, and may effectivel y radiate only 
at en e rgies in soft X-rays ba nds and above (|Sazonov et alj 
l2004| y IShemmer et al l (2006) found that a composite X-ray 
spectrum of 21 z > 4 quasars had an effective spectral index 
in the X-ray of s = 0.95, and an intrinsic absorption column 
upper limit of Nn 5. 6 x 10 22 cm " 2 , consistent with the mean 
results from I Just et alj ([2007) at lower redshift. 

Motivated by the suggestions above that the I-front 
thickness could be measurable, and by the lack of knowl- 
edge about the SEDs of high-z quasars, in this paper we 
follow the time-dependent propagation of an I-front into the 
IGM, and study its thickness for a range of possible spectra, 
quasar ages, clumping factors in the IGM, and various ion- 
izing backgrounds. Our goal is to give a rough quantitative 
assessment of the conditions under which a hard spectrum 
may be diagnosed in future o bservations. Our mode ling here 
extends the investigation of IZaroubi fc "Siikl l|2005l ). by ex- 
ploring a range of quasar spectral parameters, by improving 
the treatment of helium and of secondary ionizations by en- 
ergetic photoelectrons, and, most importantly, by following 
the time-dependent evolution of the non-equilibrium ion - 
ization structure. In an earlier study, IShapiro et al.l (|2004l ) 
used numerical simulations to study the evolution and prop- 
erties of I-fronts propagating into the high-z IGM, for three 
different source spectra appropriate for massive population 
II and III stars and quasars. Our paper extends this work 
by including recombination radiation from helium, and by 
focusing on brighter sources and covering a larger range of 
spectra (incl uding obscured spectra much harder than those 
examined bv lShapiro et all ), for which the thickness may be 
directly measurable. 

Aft er we completed this pape r, we became aware of the 
work of iThomas fc Zaroubil (|2007i ) who have recently made 
calculations similar to some of those presented here. We will 
compare our results further in §[S] 

The rest of this paper is organized as follows. In § [21 we 
describe our modeling of the I-fronts, including treatments 
of the non-equilibrium evolution, and of helium. In § [3l we 
describe the salient features - evolution, shape and size - 
of the ionization front in our fiducial model for an obscured 
source with a hard spectrum. In § [4j we present our main 
results, discussing measurements of the thickness of the I— 
Front. In § O we reiterate the inherent limitations of our 
modeling approach, before finally, in §[6] we summarize our 
results and conclusions. 



2 MODELING AND METHODS 

In this section, we describe our model for computing the 
propagation of the I-front into a uniform medium. The mod- 
eling is rather standard for the most part, but we include 
a description of our own numerical implementation here for 
completeness, and in order to describe a few details about 
helium ionization that have typically not been included in 
previous works studying cosmological quasar H II regions. 
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2.1 Ionization and Recombination 



lOsterbrock fc Ferlandl (2006) give an excellent description of 
one-dimensional spherical ionization structure calculations. 
For the most part we follow their methods, but we addi- 
tionally include time evolution (non-equilibrium ionization 
states), secondary ionization by photoelectrons, and a full 
treatment of doubly-ionized helium. 

We use the equations of one-dimensional radiative 
transfer to find the optical depth t(u) as a function of dis- 
tance from the quasar r and frequency v. The optical depth 
contribution from each species is calculated at its ionization 
threshold: 



<?o,i rii(r') dr 



(1) 



where rn{r) is the number density of species i and <7o,i is 
the photoionization cross section at the threshold frequency 
i/o,i|j The optical depth contribution at any frequency is 
then Ti(y) — n t o fi(i / )/co,i where Oi(v) is the cross section 
of species i at frequency v. The total optical depth is the 
sum of contributions from all three species 



r(v) = THl(f) + THel(^) + THellM 



(2) 



We ignore the effects of li ght travel time in cal culat- 
ing the ionization structure. As lCen fc Haimanl (|2000T ) have 
pointed out, this ends up predicting the structure that is 
observed along the line of sight. Because of the finite speed 
of light, the evolution of the gas at a distance r from the 
quasar is retarded by the interval r/c. But an observer far- 
ther along the line of sight will also receive photons from that 
point earlier by an amount r/c. So in order to correctly pre- 
dict the profile as it would be observed along the line of sight 
(for instance, in an absorption spectrum), we sh ould add the 
two e ffects, which simply cancel each other out (|White et al.l 
2003). If we wanted to predict the transverse structure as it 
would be observed, we could simply ad d in the appropriate 
delays after the evolution is calculated l|Wvithe et al.l [20051 ; 
lYull2005l ; [Shapiro et al J 1 2006b . 

For simplicity, we also ignore Hubble expansion. Even 
at z = 6, a quasar lifetime of 10 s years is much shorter than 
the Hubble time. We do note, however, that the changing 
density of the IGM could make a discernible difference to 
the shape of the ionization profile, especially for the largest 
or thickest H 11 regions, which extend to several percent of 
the Hubble radius. 

The photoionization rate of species i (per atom) is 



,-i~0,r) 

4irr 2 



(3) 



where L v is the specific luminosity, and r is the distance 
from the source. In our numerical ionization code, we tabu- 
late T; as a function of thi.o, thcI.o, and th c ii,o before cal- 
culating the evolution of the ionization structure. Added to 
that rate for hydrogen is a position-independent background 

1 Threshold ionizati on cross sections were taken from 
Osterbrock fc Ferlandl d2006h for H I and He II and from 



HTmimer*'fe***Sto"revl ll 19981") for He I. At other frequencies, we used 



the fitting formulae f r om [ Yan et al ] jl998h for He I and from 
lOsterbrock fc Ferlandl J2006h for H-like ions. 



photoionization rate, parametrized by the resulting equilib- 
rium ionized fraction Xhii. bg' 



Y 2 

1 — Xmi, bg 



(4) 



where as is the so-called "case B" recombination rate (the 
recombination coefficient to all excited levels of hydrogen, 
evaluated at T gas = 10 4 K unless otherwise specified) and 
C is the clumping factor C = (n|) / '{nn} 2 ■ This background 
arises from preexisting galaxies, and is expected to consist 
of a patchwork of H 11 bubbles surrounded by a relatively 
neutral IGM. Note that even if the ionization outside the 
quasar's H n region is patchy (with a swiss-cheese topol- 
ogy), as it is expected to be, in the interior of the H 11 
region, where the low neutral fraction results in a long 
m.f.p., it will be much more uniform (and will also have 
a higher amplitude). There will, however, be a radial pro- 
file to the backgrou nd flux, due to the clu s tering of galaxies 
around the quasar llWvithe fc Loebl 120071 : iLidz et all 120071 : 
lAlvarez fc AbelH2007f T We will first ignore these effects, as- 
suming that the I-front is located sufficiently far away that 
this bias is small, and that we are averaging the I-front over 
many small galaxy-bubbles. However, in § 14.2.41 we will re- 
turn to this issue, where we will also explore various models 
for a non-uniform ionizing background. 

In addition to the quasar's radiation and the UV back- 
ground, several other processes occur in the gas that affect 
the ionization balance and couple the ionization states of 
hydrogen and helium. 

Photoelectrons produced during ionization carry 
residual energy that can cause further ionizations 
jShull fc van Steenberd Il985h . In our calculations we in- 
clude the secondary ionization of hydrogen by photoelec- 
trons from both hydrogen and helium, and treat it as an 
on-the-spot process (the latter assumption is justified sim- 
ilarly to the argument for case-B recombination, since the 
collisional ionization cross-section i s typically larger than 
the p hoto-ionization cross-section; IShull fc van Stccnbcrg 
1985). We use the fitting formula provided bv lDiikstra et al. 
(2004) in the high-energy limit to calculate the fraction 
of energy 4>(X e ) each photoelectron will expend in further 
ionizations of hydrogen. Using the high-energy limit intro- 
duces only a small error since low-energy photons cause 
few secondary ionizations in any case and the function 
quickly approaches its asymptotic value as the photon en- 
ergy increases. Note that in the original calculations of 
<t>(x), IShull fc van Steenberd (|l985h defined x = nmi/riH 
and assumed nue/nu = 0.1 (versus 0.079 in our model), 
and XHeii = Xhii (where X Hc ii = nneii/nne)- While nei- 
ther of those conditions holds exactly in our models, they 
are reasonable approximations. To account approximately 
for the extra electrons introduced by doubly-ionized helium 
(which was not included originally) we use x = Xhii + 
("hc/^h)X Hc iii (where X H ein = nnein/nue)- It is a small 
effect in any case, since Helll is formed only where hydro- 
gen is already highly ionized, and <f> vanishes as x approaches 
unity. 

We neglect secondary ionizations of helium which hap- 
pen at less than 20% of the rate of hydrogen. The mean 
number of secondary ionizations of hydrogen per photoelec- 
tron produced in the ionization of species i is 
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n 2 = 4>{X e ) 



(Ej) - hi%i 
hv ,m 



(5) 



where (Ei) is the mean energy of photons locally ionizing 
species i: 



(Ei) 



1 



dv (hv) Oi(y) 



hv 4nr 2 



(6) 



Similarly to Ti, we tabulate (Ei) as a function of thi,o, THei.o, 
and THeii,o, before integrating the ionization structure. 

Photons produced during recombination of hydrogen 
and helium are also included in our code in the on-the-spot 
approximation. High-frequency recombination radiation is 
capable of ionizing neutral hydrogen and helium, as well as 
singly-ionized helium, so the fraction of recombination pho- 
tons ultimately ionizing each species must be calculated. In 
general the fraction of photons of frequency v that ionize 
species i is 



i(v): 



am(v) nm + cmei(y) n He i + cr H eii(^) riHeii 



(7) 



The overall fraction for a given recombination process is an 
average of yi(y) over the spectrum of the emitted radiation. 
In our calculations, we make the simplifying assumptions 
that all H I, He I and He II Lyman continua, and He II 
Balmer continuum photons are produced with the minimum 
energy. The first and last of these occur below the threshold 
for all species but HI. The He n Lyman a line is treated 
similarly, assuming all photons have the mean energy of that 
line. 

Certain higher He I lines produces in the recombina- 
tion cascade are energetic enough to ionize hydrogen, with 
an average of 0.96 H-ionizing photons ultimately produced 
for each recombination to any excit ed level of Hel (in the 
low-d ensity limit which applies here, lOsterbrock fc Ferlandl 
2006). Another source of ionizing photons is the transition 
from the first excited level to the ground state of Hell: in- 
stead of producing a He II Lyman a photon, it can instead 
occur via a two-photon process. We calculate the fraction of 
recombinations to all excited states of Hell that end in a two- 
photon emission (as a function of temperature) by interpo - 
lating the tabulated values from lHummer fc Seatonl {1964). 
An average of 1.425 photons capable of ionizing hydrogen 
and 0.737 capable of ionizin g neutral helium are prod uced 
for each two- phot on event |Flower fc Perinottd Il980l ). the 
latter of which are assumed to have the threshold ionization 
energy for helium. 

We use recom bination rate coefficients from 
iHui fc Gnedinl (|l997l ). except for the coefficient for re- 
combinations directly to the n = 2 lev el of Hell, which 
is taken from iFerguson fc Ferlandl (Il997l ). All calculations 
presented in this paper assume an isothermal gas of 
constant mean density and clumping factor. 



2.2 The Ionization Structure 

The ionization structure surrounding the quasar is calcu- 
lated by integrating the net ionization rates over time. Be- 
fore the first iteration, an array of (Ami, AhcII, AhoIii) val- 
ues, corresponding to a sequence of r values, is initialized 
to (Ahii.bg, 0, 0). The ionization and recombination rates 
are calculated as described above and used to find dXmi/dt, 



dXucii/dt, and dXneiii/dt. The integration time step is cal 
culated using 

F di 



At 



max(rHiAHi, ThoiAhcI, FhoIiAhoIi) 



(8) 



where Fat = 0.004 is just a numerical parameter optimized 
by trial and error. Decreasing Fdt by a factor of two (im- 
proving the time resolution) has no important effect on the 
results. The time step is only updated every 10 iterations, 
and to help damp numerical oscillations we actually use 
Af ncw = 0.75 At id + 0.25 At as the time step. 

In order to increase the speed and stability of the inte- 
gration, the ionized fractions for the inner part of the array 
(at low r) are frozen and no longer recalculated once a stable 
equilibrium is achieved. Every tenth iteration we calculate 
the maximum radius at which the ionization states are in 
equilibrium, defined as r eq = max(r) for which 



Am 



> 5 x 10 7 At 



V dt 

X Hc „(^i)">5xl0 7 A i 

THI.O < 1000 



(9) 

(10) 
(11) 



The last condition on the optical depth is introduced to ex- 
clude regions far ahead of the front, where time scales grow 
long because the radiation field is so weak. For runs with a 
high uniform ionizing background we decreased the optical 
depth limit to thi.o < 100 because of the reduced opacity. 
Once r eq is found, the next iteration is begun, and the ion- 
ized fractions at r > r cq are updated (and constrained to 
the interval [0, 1]). Using our fiducial parameters (described 
below and in § [3j , it took 76000 iterations to evolve the 
structure through 5 x 10 s years. 



2.3 Numerical Parameters 

As mentioned above, ionization rates and mean photon en- 
ergies are tabulated as a function of the threshold optical 
depths of each species. The quantities are evaluated on a 
grid of logarithmically spaced values of thi.o, i~Hei,o, and 
T"Heii,o- The dimensions of the ionization rate grids are 332 
by 301 by 301 values, while the mean photon energy grids 
are 416 by 61 by 61 values. Each optical depth range ex- 
tends from 10 -8 up to the maximum optical depth for which 
Ti/ max(rj) > 10~ 6 , which depends on the source spectrum. 

The simulations presented here use a grid of 1300 
equally spaced radius values, extending from 5 x 10 4 pc (all 
distances in this paper are in proper, not comoving, units 
unless stated otherwise) to 1.5 -Rstrom, where iistrom is the 
radius of the classic equilibrium Stromgren sphere, 



iistr 



(12) 



Note that this expression ignores helium and secondary ion- 



= 2 x 10" 



and C = 1, i?s 



trom — 



izations. With 
22.2 Mpc at z = 6. 

We tested for convergence by doubling the resolution, 
which had no effect on the results. For our canonical param- 
eters the step size works out to Ar = 2.56 x 10 -2 Mpc with 
an outer radius of 33.4 Mpc. 
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We have tested our code several ways to ensure reason- 
able numerical performance. As mentioned above, we tested 
for convergence by doubling either the time or spatial res- 
olution, neither or which had any important effect on the 
results. Using a soft spectrum (so that the front is sharp) 
we find very good agreement with analytical estimates for 
ionized region size and front velocity. We also find agreement 
with the Cosmological Radiative Transfe r Codes Compari- 
son Project ( "test 1" in llliev et al"1l2006l) calculations, and 
the examples presented in lOsterbrock fc Ferlandl (|2006h . 



2.4 Physical Parameters 

The temperature of the gas is important, because it controls 
the recombination rates. All calculations presented in this 
paper assume an isothermal gas. We justify this assumption 
by noting that most of the photoionization heating occurs 
somewhat ahead of the front where the ionized fractions are 
small and recombination rates are therefore low, whereas 
the tempe rature inside the front is f airly uniform (a s can 
be seen in iBolton fc Haehneldl2007bl ; llliev et alj|2006h . We 
experimented with introducing artificially varying tempera- 
ture profiles in the outer regions of the front. We made the 
temperature a function of optical depth so that the tempera- 
ture variations would be tied to the location of the front. We 
tried several profiles, from smooth func tions mimicking the 
temper ature profiles in the appendix of IBolton fc Haehneltl 
(2007b), to a step function increasing the temperature by 
a factor of 5 beyond thi.o > 0.001. We found that none of 
these had a significant effect on the profile, as long as the 
variation occurred outside of the region that had reached 
equilibrium. The ionization rate totally dominates the re- 
combination rate in the moving part of the front, reaching 
parity only as the ionized fraction approaches equilibrium. 
We conclude that since the heating of the gas occurs before 
it reaches ionization equilibrium, and since the exact recom- 
bination rate is irrelevant outside of the equilibrium region, 
changes in the recombination rate due to temperature vari- 
ations will not be an important factor in determining the 
shapes of the fronts 

We assume the number density of hydrogen atoms sur- 
rounding our sources is the mean density of the IGM at 



z = 6, n H = 2.2 x 10~ 7 cm~ 3 



(1 + *) 



7.6 x 10" D cm" 



The 



helium mass fraction is 0.24, yielding nn c = 6.0 x 10 cm . 
The electron temperature is T gaa = 10 4 K. For our fiducial 
case we use a clumping factor of C — 1. 

The spectra we have used are absorbed power laws, as 
show in Figure [T] with specific luminosity 



exp i 



<ym\y) H OYL&yv) 



(13) 



where s is the spectral index and Ah and iVke are the column 
densities of hydrogen and helium intrinsic to the source. The 
spectrum is normalized to produce ionizing photons at rate 
Q before the absorption, so the ionizing photon luminosity 
after absorption Q a bs can be much lower. Note that high 

2 Changes in temperature due to radiative transfer can also have 
dynamical effects on the gas around the I-front; but in the limit 
where the recombination rate is unimportant, our conclusions on 
the ionized fraction should also be insensitive to such dynamical 
effects. 




1000 



E eV 



Figure 1. An absorbed power law spectrum (arbitrary units) 
with s = 1.6 and TVjj = 10 19,2 cm -2 (solid curves) plotted in 
linear (upper panel) and logarithmic (lower panel) flux units. The 
dashed lines show the unabsorbed spectrum. The vertical dot- 
dashes lines indicate the hydrogen ionization threshold at 13.6eV. 



energy photons can cause multiple secondary ionizations, so 
Qabs is not exactly the total ionization rate the spectrum 
produces. It is difficult to take secondary ionizations into 
account a priori, since (as we discussed earlier) the number 
of secondary ionizations is a strong function of the local free 
electron density, and therefore depends on where the photon 
is absorbed. Other factors can have an even larger effect on 
the ionizing efficiency of the spectrum, such as the fact that 
the hardest photons can escape the region entirely. 

As mentioned in the Introduction, the ionizing spectra 
of quasars at z > 6 is uncertain. For our fiducial model for 
an obscured source, we have chosen Q = 2 x 10 57 s _1 (corre- 
sponding roughly to the expected value for th e z ~ 6 quasars 
in the SDSS; e.g. iMesinger fc Haimanlliooj ). s = 1.6, and 
log^H/cm -2 ) = 19.2 (which results in a post-absorption 
photon luminosity of Qabs = 8.3 x 10 55 s _1 ). We will con- 
sider variations in the power-law index and obscuring col- 
umn density that are consistent with expectations based on 
the lower redshift data. 



3 EVOLUTION, SHAPE AND SIZE OF THE 
IONIZATION FRONT 

In this section, we present our results for the time-evolving 
I-front around a fiducial source with hard spectrum, and 
discuss the broad features that can be relevant to future 
measurements. Quantitative results will be presented and 
discussed in the next section. 



3.1 The Time-Evolving Ionization Front 

An often used proxy for the thickness of an ionization 
front around a quasar is the mean free path Zmfp = 
[ nm cm ((f)) ]~\ where (y) is a characteristic frequency 
of the ionizing photons. The thickness can be estimated, 
for example, by adopting nm = 0.5 tih, and the mean fre- 
quency of ionizing photons emitted by the quasar for (is). 
The resulting value of Zmfp, however, is an accurate proxy 
for the thickness of the I-front only for particular definitions 
of the thickness, and only under particular conditions. For 
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instance, the ionizing spectrum is hardened as the photons 
travel outward through the gas, since low-energy photons 
are more readily absorbed. As a result, the outer edge of the 
front will be spread over a longer path than the inner edge. 

Another complication is that the ionization fronts 
around short-lived quasars are expected to be propagat- 
ing outward (rather than corresponding to the static edge 
of an equilibrium Stromgen sphere), which has a strong 
effect on the shape of the front. It is widely noted (e.g., 
IShapiro fc Girouxl Il987n that the ionization state of the 
gas surrounding a typical high-redshift quasar will not have 
time to reach equilibrium with the radiation field during a 
quasar lifetime of 10 7 to 10 s years0 since the time scale 
to reach equilibrium is roughly the recombination time 
t rec = (Cqshh)" 1 = 1.62 x K^C -1 [(l + z)/7]~ 3 years. 
This timescale can also be thought of as the time required for 
the source to emit one photon for each hydrogen atom within 
the Stromgren sphere. Throughout the process of establish- 
ing equilibrium, the shape of the ionization front is changing, 
therefore we must take into account non-equilibrium effects 
in order to correctly predict the ionization front thickness. 

Figure [2] shows the evolution of the ionization struc- 
ture around our fiducial absorbed power-law source. The 
hydrogen front starts out quite thin, but becomes thicker 
and thicker as it propagates outward. When the source is 
turned on, the gas closest to it is very quickly ionized up to 
its equilibrium level, while gas a little further away is ionized 
much more slowly, due both to the intervening absorption 
and geometric dilution of the radiation. This means that the 
inner edge of the front is evolving faster than the outer edge, 
which effectively tilts the Xhi curve, making it steeper. This 
creates a very thin transition region at first, which thickens 
as the front propagates outward and slows down. Depend- 
ing on how the thickness of the front is measured, it can 
start out thinner than the mean free path, and end up much 
thicker. 

These effects can be seen even with a monochromatic 
spectrum, but, as illustrated in Figure [3] a spectrum with 
a broad continuum of ionizing radiation can enhance the ef- 
fect. Early on (e.g. at 10 6 years), photons with lower energy 
(and therefore high ionization cross section) will dominate 
the ionization at the front, while higher-energy photons es- 
cape to larger distances and pre-ionize the medium (as seen 
in the long tail of curves representing the ionized fraction). 
This makes the front relatively thin because of the short 
mean free path of these photons. As the front moves outward 
(e.g at 10 s years), the lowest energy photons will continue 
to be preferentially absorbed closer to the source (and will 
dominate the ionization balance in the inner regions), leav- 
ing the higher-energy photons to dominate in the outer parts 
of front. This hardened spectrum creates a thicker front be- 
cause of the longer mean free path. 

Based on the above features, it will be useful for the dis- 
cussion below to divide a snapshot of the ionization structure 
at a time t < t lac into three distinct regions: Closest to the 
quasar, in the equilibrium zone, the ionization rate has al- 



3 A combination of overdensity and high clumping factors in the 
typical quasar environment could, however, reduce the recombi- 
nation t ime to a value perhaps as short as ~ 5 X 10 6 years at 
z = 6.4 <Yu fc Lull20"u5h . 
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Figure 2. Ionized (dashed red) and neutral (solid blue) hydrogen 
fractions (a); and neutral (solid cyan), singly- ionized (dot-dashed 
magenta) and doubly-ionized (dotted yellow) helium fractions (b) 
versus radius (proper coordinates). Curves in the upper are for 
quasar ages of 4 < log(t/years) < 10 in 1 dex increments. In the 
lower panel, only log(t/years) = 6, 8, and 10 are shown for clarity. 
By 10 10 years the entire ionization structure is in equilibrium. 
Model parameters are Q = 2xl0 57 s _1 ; s = 1.6; log(7VH/cm -2 ) = 
19.2. The mean photon energy for this spectrum is s=s 240 eV, and 
the corresponding mean free path is ~ 3 (proper) Mpc. 



ready reached equilibrium with the recombination rate. Just 
outside the equilibrium region is a rapid-ionization zone, 
where relatively neutral gas is suddenly being exposed to in- 
tense ionizing radiation. Further out there is a pre-ionization 
zone, where the relatively rare high-energy photons pene- 
trate and begin to ionize the gas, but recombination rates 
are low because of the low ionized fraction. 

Figure [4] shows a snapshot of the local mean ionizing 
photon energy (equation |6}, and the corresponding mean 
free path (calculated using the local neutral fraction), ver- 
sus neutral hydrogen fraction, for the profile at 10 s years 
shown in Figure [2] The spectrum is relatively unaffected by 
hardening up to Xhi ~ 0.1, but the mean energy increases 
rapidly after the neutral fraction reaches 0.3. Similarly, the 
mean free path at first decreases in inverse proportion to 
the increasing neutral fraction, but reaches a minimum at 
Xhi ~ 0.3 and begins to grow as the spectrum hardens. This 
explains why the outer tail of ionization is so thick compared 
to the inner face of the front. If we measure the thickness of 
the front from Xm = 0.1 to 0.9, we find (as shown in the fig- 
ure) that it is barely consistent with the range of Imfp over 
the same interval: both approaches to estimating the Zmfp 
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r (Mpc) 

Figure 3. Ionization structure for an absorbed power law spec- 
trum (solid, same parameters as Figure [2]| and for a monochro- 
matic spectrum (dashed). Ionized (red, decreases with r) and 
neutral (blue, increases with r) hydrogen fractions are shown at 
t = 10 6 years (left set of curves) and t = 10 8 years (right set of 
curves). The monochromatic spectrum is normalized to produce 
ionizing photons at the same rate Q a b s = 8.3 X 10 55 s — 1 and 
has a photon energy of 89 eV, which is the local mean ionizing 
photon energy (weighted by cross section) (-Eh) of the absorbed 
power law spectrum. Note that a convenient way to estimate the 
thickness of the front is to compare where the ionized and neutral 
fraction curves cross log(X) as —1 at a given time. 




log(X HI ) 

Figure 4. The upper panel shows the mean ionizing photon en- 
ergy (En) for hydrogen (solid curve), neutral helium (dashed 
curve), and singly- ionized helium (dot-dashed curve). The lower 
panel shows the local mean free path versus neutral hydrogen 
fraction at 10 s years (solid curve; using the same model as in 
Figure(2j. The lower dashed line shows £mfp calculated using the 
r = value for the local mean ionizing photon energy (weighted 
by cross section) (-Eh) = 89 e V, while the upper dashed line is 
calculated using the r = overall mean ionizing photon energy 
(E) = 239 eV. The horizontal dot-dashed line indicates the thick- 
ness of the front as measured from Xm = 0.1 to 0.9. 



based on an average photon energy of the source become 
inaccurate because of the spectral hardening. The intersec- 
tion of the three curves at log(Xni) = —0.1 is merely a 
coincidence. Indeed, recall that the ionization structure is 
not in equilibrium, so the width is constantly changing. The 
estimate of the m.f.p. based on the mean photon energy of 
239 eV and a neutral medium yields ~ 3 Mpc, whereas the 
path corresponding to the range 0.1 < Xui < 0.9 in the 
equilibrium case is much longer, ~ 15 Mpc. 



log(t/ year) 

Figure 5. Width measurements versus time for the ionization 
front shown in Figure [2] The solid curve shows the thickness of 
the front as measured from Xhi = 0.1 to 0.9, while the dashed 
curve shows the thickness between -Xjji = 10 — 3 and 10 — 2 5 . 




r (Mpc) 

Figure 6. Ionization structure with (solid) and without (dashed) 
ionization of hydrogen by helium recombination photons. Ionized 
(red, decreases with r) and neutral (blue, increases with r) hydro- 
gen fractions are shown at t = 10 6 years (left set of curves) and 
t = 10 8 years (right set of curves). Note that this recombination 
radiation accounts for almost all of the influence of helium on the 
hydrogen ionization structure. 



Figure [5] shows two different measurements of the thick- 
ness of the front versus the logarithm of the time. The def- 
initions of these thicknesses will be justified in § 21 but for 
now we note that the outer part of the front grows thicker 
all the way up to about i roc , while the inner part of the front 
comes into equilibrium much sooner. 

The presence of helium also has important effects on 
the hydrogen ionization structure, as can be seen in Fig- 
ure [6] The influence is due to the coupling of the ionization 
states by recombination radiation. By reprocessing high- 
energy photons (which, in the absence of helium, would es- 
cape to large distances) into lower-energy recombination ra- 
diation (which is absorbed on the spot, mostly by hydrogen), 
helium lowers the hydrogen neutral fraction in the equilib- 
rium zone. The difference is a factor of ~ 3 in the region 
where Xm < 10~ 2 at 10 s years; this region is relevant in the 
analysis of the proximity zone of the quasar in Lya//3 ab- 
sorption. In particular, neglecting helium can result in pre- 
dicting a much thicker front, depending on how thickness is 
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12 3 4 



r (Mpc) 

Figure 7. Ionization structure with (solid) and without (dashed) 
secondary ionizations by high energy photoelectrons. Ionized (red, 
decreases with r) and neutral (blue, increases with r) hydrogen 
fractions are shown at t = 10 6 years (left set of curves) and t = 
10 s years (right set of curves). 




1 2 3 4 5 



r (Mpc) 

Figure 8. Ionized and neutral hydrogen fractions for a source 
with no intrinsic absorption (N-r = 0, dashed) or with 
log(iV H /cm- 2 ) = 19.2 (solid). Spectra are normalized to pro- 
duce the same number of photons Q a b s = 8.3 X 10 55 s . Other 
spectral parameters are the same as in Figure [2] 



measured (and a factor of several too large Rp/ R a ratio; see 
below). 

Secondary ionizations have a similarly important effect, 
but on the outer tail of the ionization structure, as shown 
in Figure [7] Neglecting secondary ionizations can reduce the 
thicknesses we would predict in the outer part of the front 
by several Mpc. 



3.2 The Effects of Spectral Hardness 

As evident from the above discussion, the shape and hard- 
ness of the spectrum can have dramatic effects on the ion- 
ization structure. Generally a harder spectrum will produce 
thicker fronts, although as we shall see below, this is not al- 
ways the case for certain definitions of front thickness. The 
mean photon energy is not enough to characterize the shape 
of the ionization structure it will produce. For instance, even 
a relatively weak high-energy component can produce an en- 
hanced pre-ionization region, while a low-energy component 
can increase the ionized fraction of the equilibrium region 
but have little effect on the thickness of the rapid-ionization 
zone. Enhancing the effect of a high-frequency component is 
the fact that energetic photons can cause multiple secondary 
ionizations, especially in a relatively neutral medium. 

A harder spectrum will also generally produce a "shal- 
lower" ionized region, in the sense that the ionized fraction 
will be lower in the inner part of the ionized region since the 
ionization cross section is lower for high-energy photons. At 
the quasar ages we are considering this effect can be consid- 
ered to affect mainly the equilibrium region. 

Two parameters control the shape of the spectrum in 
our models. Increasing the spectral index s makes the spec- 
trum softer, while increasing the intrinsic hydrogen column 
density A^h cuts off more low-energy emission, hardening the 
spectrum. 

For reference, Figure [8] shows the fronts produced by 
a spectrum with no intrinsic absorption (JVh = 0). As ex- 
pected, the fronts are much sharper than those in Figure 

m 



3.3 The Size of the Ionized Region 

Before discussing the thickness of the I-fronts in detail in 
the next section, it is useful to discuss the overall size of the 
H II region. As mentioned in the Introduction, the size of 
the H II regions could be probed either in 21cm studies, or 
by Lyman line absorption spectra of individual sources. The 
former technique is sensitive to Xm ~ 1, whereas the latter 
technique (with typical deep Keck spectra) probes the much 
lower neutral fraction Xhi ~ 10~ 3 . In Figure|5] we therefore 
show the size of the ionized region measured out to ro.5, 
where Xm = 0.5 (top panel), and to r_3, where Xm = 10 -3 
(bottom panel). The shaded plots show the H II region sizes 
at 3.16 x 10 7 years (corresponding roughly to the Salpeter 
time, and to the expected ages of bright quasars; see, e.g., 
lMartinill200i ). as a function of the intrinsic absorption Nn 
and spectral index s. 

Remember that the softest spectrum is in the lower 
right corner of the graph, while the hardest spectrum is in 
the upper left. The trend in both panels is for the radius 
to decrease with increasing Ah, as Q a b s decreases due to 
the absorption, and as more of the photons escape from the 
ionized region because of their longer MFP. Decreasing the 
spectral index s (which also produces a harder spectrum), 
on the other hand, results in a larger H II region, since the 
harder spectrum is less attenuated by the intrinsic absorp- 
tion, meaning Q a b s is larger. 

Considering both H II region size measurements to- 
gether leads to an interesting conclusion. Several re- 
cent papers (iMesinger fc Haiman 2004 ; Maselli et al.ll2007l ; 



iLidz et all 120071 ; iBolton fe Haehneltl l2007aT ) have explored 
the accuracy of R a and Rp, defined as the radii where the 
flux in the redshifted Lya and Ly/3 region of the spectrum 
falls below some fixed threshold, as a proxy for the loca- 
tion of the I-front. They found that R a and Rp are biased 
measurements, underestimating the "true" distance to the 
I-front (which we might choose to locate at ro.5) by typically 
20-30%. Here we discover another effect to worry about: as 
the spectra get harder (and especially as the obscuring col- 
umn increases), the radius of the highly ionized region r_3, 
decreases much more quickly than the half-ionization radius 
ro.5. This is because the harder photons are ineffective at 
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Figure 9. The size of the ionized region in Mpc out to Xhi = 0.5 
(ro.5, a) and Xjji = 0.001 (r_3, b) versus the intrinsic absorption 
Nn (in cm -2 ) and the spectral index s at 3.16 X 10 7 years and 
with Q = 2x 10 57 s — 1 . Each block represents one simulation, with 
the size indicated in numerals and by the shading. The key to the 
shading appears at the top of each graph. 

maintaining a high ionized fraction because of the low ion- 
ization cross section, instead tending to escape into more 
neutral gas farther out. This produces a thicker front and a 
more neutral H n region. For hard spectra, this significantly 
increases the bias, with the underestimate reaching values 
as high as a factor of ~ 8 for the spectra with the largest 
obscuring column. 



4 MEASURING THE THICKNESS OF THE 
I— FRONT 

4.1 Front Thickness in 21cm and in Lyman 
Absorption Spectra 

As mentioned above, we are interested in two different re- 
gions of the ionization front: the highly ionized inner face 
of the front at 10 -3 < Ahi < 1CP 2 ' 5 which is approximately 
the region probed by Lyman-series absorption spectra, and 
the less-ionized outer tail of the front at 0.1 < Ahi < 0.9, 
which is the region potentially accessible to redshifted 21 
cm observations. 



iMesinger fc Haimanl (|2004l ) proposed the idea that the 
difference between Rp and R a , as defined in the previous 
section, could be useful to probe the spectral hardness, with 
the two values incre asingly separated for harder spectra. In 
a subsequent study, iBolton fc Haehneltl (|2007al ) simulated 
a large set of quasar Lyman a and (3 absorption spectra 
and studied the possibility of using the transmission win- 
dow created by the ionized region around the quasar to con- 
strain the mean ionization fraction of the IGM. They found 
that density fluctuations induce a large scatter in their mea- 
sured quantity, which is the ratio of radii inferred from each 
spectrum R/s/R a . However, according to their simulations, 
Rft/Ra — 1.2 would be distinguishable from Rp/R a = 1 at 
the 2a level with about 20 high-quality quasar spectra. As 
a proxy for this ratio, we use r_2. sA-3, the ratio of radii 
at which Am = 10 -2 ' 5 and Ahi = 10~ 3 , which correspond 
roughly to the points where flux would be lost in the Ly-/3 
and Ly-a absorption spectra, respectively, for a typical deep 
Keck spectrum of a z — 6 quasar. 

In contrast to the Lyman line absorption spectra, the 
21cm observations will be sensitive, at z — 6— 10, to features 
corresponding to neutral fractions of order unity, with re- 
gions that have Ahi J$ 0.1 corresponding to "holes" in tomo- 
graphic 21cm maps. Although the effective spatial resolution 
of such maps is highly dependent on the specific interfero- 
metric instrument and its configuration, a rough target value 
for next -generation facilities may be ~ 1 p hysical Mpc (see, 
e.g., the review by iFurlanetto et alj 1200a ). We adopt here 
do.u the difference between the radius at which Ami = 0.1 
and Ahi = 0.1, as a proxy for the thickness of the I-front 
that may be measurable with 21 cm instruments. 

Figure \W\ depicts our measures of the ionization front 
thickness around a quasar with an age of ~ 3.16 x 10 7 years, 
on the same grid of parameters used in Figure [9] The first 
panel shows do.i, which exceeds our estimated detection 
threshold of ~ 1 Mpc for a large region of the parameter 
space. The thickness increases with increasing spectral hard- 
ness up till log(AH/cm -2 ) = 19.8, after which is decreases 
as Ah increases, due to the overall shrinking of the ionized 
region discussed in § 13.31 and shown in Figure [9] This will 
be the pattern for most of the grids we will discuss: the do.i 
widths will always increase with decreasing spectral index, 
and will grow with increasing absorption up to a point, but 
will eventually shrink for large enough Ah. 

The r_2.5/r_3 ratio in the second panel is almost to- 
tally insensitive to the spectral index. Making the spectral 
index smaller does make the front slightly thicker, but it 
also moves the front outward by a small amount, which keeps 
the ratio constant. The ratio increases dramatically between 
log(A H /cm -2 ) = 18.6 and log(A H /cm~ 2 ) = 19.2, becoming 
large enough to be potentially detectable (according to th e 
criterion defined above based on lBolton fc Haehneldl2007al ) ■ 
thanks to the fact that for the harder spectra r_3 reaches 
equilibrium faster, while r-2.5 continues to move outward 
a little longer, which stretches out the ratio. The intrinsic 
absorption must be at least log(AH/cm -2 ) = 19.2 for the 
ratio to exceed 1.1. The ratio exceeds 1.2 for any spectral 
index at log(A H /cm~ 2 ) > 19.8. 

Note that the contours for the two size measures in the 
Ah-s plane have different orientations. This is not surpris- 
ing, since variations of s and .Ah can affect the profile differ- 
ently at different locations in the front. This suggests that a 
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Figure 10. Ionization front thickness measures at 3.16 xlO 7 years 
as measured by do.i (in Mpc, top) and r_2.s/r_3 (bottom). 



Figure 11. Widths at 1 X 10 8 years as measured by do.i (in Mpc, 
top) and r_2.5/r_3 (bottom). 



combination of 21cm and Lyman-series absorption measure- 
ments could break the degeneracy inherent in the separate 
observations and allow both parameters to be constrained 
simultaneously. 



4.2 Parameter Variations 

We next examine how our results change when we vary other 
model parameters. 



4-2.1 Source Lifetime 

In Figure 1111 we show the same thickness measures as in 
Figure 1101 but for quasar three times older, at t — 10 s 
years. The figure reveals that the outer tail of the front 
(do.i) has generally gotten thicker, but otherwise follows 
the same trends we saw at 3.16 x 10 7 years. The inner face 
ratio r_2.5/r_3 has remained small and almost unchanged 
for the softest spectra with JVh < 10 18 ' 6 cm -2 . These spec- 
tra create steep fronts and deep equilibrium regions with 
Xhi < 0.001, so even at 10 8 years r-2.5 /rs is still measur- 
ing the steep slope of the rapid ionization zone. There is a 
sharp increase in the front thickness at all spectral indices 
between JV H = 10 186 and 10 192 cm~ 2 . This is because for the 
harder spectra the equilibrium region reaches Xhi > 10 -3 



(but is usually still below Xhi = 10 -2 ' 5 ) before 10 8 years, 
meaning that r_3 is the distance to a fixed point in the 
equilibrium zone, while r-2.5 is the distance to a point in 
the still-receding ionization front. This produces the high 
ratios seen in the figure. 



4-2.2 Source Luminosity 

Reducing the (pre-intrinsic-absorption) ionizing photon lu- 
minosity from Q = 2 x 10 57 s" 1 to 2 x 10 56 s" 1 produces the 
widths shown in Figure 1121 Reducing the luminosity by a 
factor of 10 reduces the time scale by a factor of 10 and the 
Stromgren radius by a factor of lO 1 ^ 3 = 2.15. The shape of 
the pre-ionization and rapid-ionization zones (since they are 
largely unaffected by recombinations) are therefore almost 
identical to the Q = 2 x 10 57 s^ 1 models at 3.16 x 10 6 years 
(not shown), when the same total number of photons has 
been emitted. And indeed the do.i values are very close, and 
consequently are lower than the Q = 2 x 10 57 s _1 models at 
3.16 x 10 7 years by a factor of about 1.5 to 2. The same rea- 
soning does not apply to the r-2.5 /f-z ratio, which end up 
with similar values. The lower luminosity means that r_3 is 
frozen into the equilibrium zone earlier, but r_2 also moves 
out more slowly, so the ratio ends up similar to the higher 
luminosity case. 
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Figure 12. Widths with Q = 2 X 10 56 s _1 at 3 X 10 7 years as 
measured by do.i (in Mpc, top) and r-2.5 /V— 3 (bottom). 



Figure 14. Widths with C = 5 at 3.16 X 10 7 years as measured 
by do. 1 (in Mpc, top) and r_2.5/r_3 (bottom). 
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Figure 13. Ionized (red) and neutral (blue) hydrogen fractions 
for C = 1,5,20 (bottom-to-top for the neutral fraction) at t = 
10 s years. Spectral parameters are the same as in Figure [2] 



ure 1131 Increasing the clumping factor also decreases the 
recombination time and Stromgren radius. With C = 5, 
tree = 3.2 x 10 s years, while with C = 20, t ICC = 8.1 x 10 7 years 
(both at z = 6), meaning that the ionization structure could 
approach equilibrium within a quasar lifetime. 

Figure [14] reproduces the grids of spectral parameters 
with a clumping factor of C = 5. The thickness measures 
follow the same trends seen with (7=1. Increasing the 
clumping factor produces almost the same grid of do.i values 
because it has little effect on the high-energy photons that 
pre-ionized the outer region. 

The r_2.5/V_3 ratios on the C = 5 grid are quite a 
bit higher than with (7 = 1. The ratio exceeds 1.2 for all 
spectra with Nu > 10 18 ' 6 cm -2 . A clumping factor of (7 = 
20 (also not shown) produces a very similar grid for the 
more obscured spectra, and produces even higher ratios for 
log(AT H /cm- 2 ) < 19.2. 



4-2.3 Clumping Factor 

The clumping factor effectively increases the recombina- 
tion coefficient, which has the effect of raising the neutral 
fraction in the equilibrium zone, but has little effect on 
the rapi d- ionization or pre-io nization regions, (as empha- 
sized in ICen fc HaimanlboOOr ), and as can be seen in Fig- 



4-2.4 The Ionizing Background 

The IGM is expected to have already undergone some non- 
negligible ionization by galaxies, and possibly by fainter 
quasars, by the time the massive black holes producing 
the observable z > 6 quasars appear. The ionizing back- 
ground may have both a smooth X-ray component (e.g. 
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Figure 15. Widths with a uniform ionizing background, 
-^hii, BG = 0.75. Note the here we are using the distance from 
-Xhii = 0.9 to Xjjii = 0.8 to measure the outer front (top) and 
still using r_2.5/r_3 (bottom) for the inner front. Both are at 
3.16 x 10 7 years. 



IVenkatesan et alj|200ll ; IOh]|200ll) from pre-existing fainter 
quasars, and a "swiss-cheese"-like component of smaller ion- 
ized bubbles surrounding (clusters of) galaxies, the latter of 
which may be highly clustered aro und the massive halos 
hosting our quasars of interest (e.g. | Wvithe fc Loebl 120071 : 



iLidz et all [20071 ; I Alvarez fc Abell I2007b 7 

Here we first calculate a grid of thickness measures with 
a purely uniform background Tbg = 4.4 x 10~ 17 s _1 (see the 
discussion of radiation transfer below), producing an equi- 
librium Ami, bg = 0.75. We will then study three variations 
of this model, to address the following effects: (i) modifica- 
tion of the background flux due to radiative transfer effects 
caused by the ionization by the quasar, (ii) a non-uniform 
(swiss-cheese) background ionization topology, and (iii) en- 
hanced ionization near the quasar due to clustering of pre- 
existing galaxies. 

In the presence of a background, there are fewer neutral 
hydrogen atoms to ionize, and the ionization fronts travel 
farther in a given amount of time. The lower neutral hydro- 
gen density also means the m.f.p. is longer, so the fronts are 
somewhat thicker. The equilibrium neutral fraction is rel- 
atively unaffected as long as it is well below 1 — Ami, .bg- 
Setting Xhii, bg = 0.75 obviously makes our previously- 



defined measure for the thickness, do.i, unusable, so we here 
instead examine the distance from Ami = 0.9 to Ami = 0.8 
(and we can still use r_2.s/r_3). Both of these thickness- 
measures are shown in Figure [T5l again at 3.16 x 10 7 years, 
as a function of spectral slope and obscuring column density. 

The new measure of the outer part of the front behaves 
much like do.i did for the grids with no background (ex- 
cept that it is obviously smaller overall). The figure also 
shows that r_2.s/r_3 is somewhat larger than for the no- 
background case. The reduced neutral hydrogen density al- 
lows the front to move more quickly, but has little affect on 
the equilibrium ionized fraction where Am C 1 - Ami, .bg- 
This means that the inner (most highly ionized) edge of the 
front slows as it approaches the equilibrium value, while far- 
ther out the front is still receding at a faster rate, effectively 
thickening the front. 

There is an inconsistency in the way the background 
ionization rate was added to the quasar flux above, however 
(which is shared by many other similar studies in the litera- 
ture). We specified an ionization rate per neutral hydrogen 
atom r_BG which is constant, equivalent to a uniform back- 
ground flux that does not depend on the ionization state of 
the gas. In fact, the flux from a homogeneous population 
of emitters will depend on the mean free path, and hence 
on the ionization state of the gas. The more highly ionized 
the gas is, the farther ionizing photons will travel through 
it, meaning that more distant sources can contribute to the 
local flux. In fact, the mean free path within a highly ion- 
ized H n region might exceed the size of the region, meaning 
that the background will be the sum of the flux from all the 
sources inside the region, and the background will grow as 
the region expands. 

In a homogeneous medium with a total photon emission 
coefficient of j (photons per unit time per unit volume), the 
total incident flux (integrated over all angles) is simpljjj 



J = 



Tin (I - Xnil)(T 



= j h 



(14) 



where a is the absorption cross section for the photons. Since 
r = Jo, the background ionization rate should clearly de- 
pend on the ionization state of the gas. Similarly, the flux 
in the center of a sphere of radius r (where ?mfp r, and 
j = outside the sphere) is J = j r. 

In our model the medium is obviously not homogeneous, 
since Ami varies with position and time (and later we will be 
interested in a variable j as well) . Solving for the background 
flux J exactly in an inhomogeneous medium requires com- 
putationally expensive 3-D radiation transfer (and still 2-D 
in the azimuthally symmetric case), so we have implemented 
an approximate radiative transfer algorithm to estimate the 
effect using a 1-D calculation. We assume all of the back- 
ground photons are at 13.6 eV and calculate the flux using 
the equation: 



J(r) = J (r) + Ji„(r) + J out (r) 

J (r) = j(r) Wp(r) [1 - exp(— Ar/luw(r))] 



(15) 
(16) 



4 Usually an emission coefficient is defined as the amount of 
power emitted per unit solid angle per unit volume. Here we use 
a somewhat different definition which is convenient for the task 
at hand. 
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Figure 16. Comparison of neutral fraction with no ionizing back- 
ground (solid), with a uniform background without background 
radiation transfer (dashed), and the same background with ap- 
proximate RT (dotted) at log(t/years) = 6, 7, and 8 (left to right). 
Spectral parameters are the fiducial ones used in Figure[2]and the 
backgrounds are chosen to produce Xhii, BG = 0.75. Note that 
the RT effects make a noticeable difference in the 10 8 year case: 
the reduction of the mean free path by the quasar boosts the 
contribution of background galaxies to the total ionization rate 
inside the H II region. 

J in (r) = ^J^ j(r') (^j eNC-O-roW]^/ (17) 
1 f°° i 

J out (r) = - j( r ') e I-oW-ro(r )] dr , (lg) 

^ Jr+Ar 

where the first term is the flux from within one radial step 
Ar, the second term is the estimated flux from within r — 
Ar, the third term is the estimate flux from outside of r + 
Ar, and ro(r) is the optical depth at 13.6 eV as defined 

by equation [T] The ad hoc factor (^t) * s use< A to decrease 
the contribution to the flux from spherical shells at small r 
(roughly approximating the decrease in the total emitting 
volume in 3D at small radii). At r ^> Zmfp with uniform 
j and -Xjni, equation [T5] reproduces the analytical result 
(equation I14|l . 

Figure \W\ compares the ionization structures with and 
without the approximate RT algorithm. The photon emis- 
sion coefficient is j = as "h^hubg = 8.3 x 10 -22 s -1 cm -3 . 
For quasar ages of < 10 7 years, the radiation transfer makes 
no difference to the ionization structure. The central ionized 
region is simply too small to allow a large additional back- 
ground to build up. At 10 8 years, RT decreases the neutral 
fraction, but only for —4 < log(XHi) < — 1.5. This is because 
the ionized region has grown large enough for the diffuse 
background to become significant, but the quasar flux still 
dominate close the quasar, and the MFP is still very small 
for 13.6 eV photons once the neutral fraction gets large, 
so the extra flux is confined to the inner part of the front. 
We conclude that radiation transfer of flux from a uniform 
ionizing background could be important for Lyman-series 
observations of old quasars or with high backgrounds, but 
will be less important for 21 cm observations. 

Next, let us examine the effect of a patchwork of pre- 
existing ionized bubbles along our line of sight. This is the 
type of background profile that results from the "swiss- 
cheese" topology of small ionized bubbles (generated by 
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Figure 17. Comparison of the neutral fraction with no ionizing 
background (solid), with a smooth ionizing background (and ap- 
proximate radiation transfer) chosen to produce Xhii, bg = 0.75 
(dashed) and with a "swiss-cheese" background (dotted) of 1 Mpc 
ionized bubbles separated by 0.79 Mpc neutral regions, which 
would produce a similar average ionization. Spectral parameters 
are the fiducial ones used in Figure [2] The panels show the ion- 
ization structures at 10 7 , 3 X 10 7 , and 10 8 years, respectively. 



neighboring galaxies and possibly AGN) embedded in a 
neutral IGM. Naively, we might choose to fill 75% of the 
IGM with ionized bubbles in order to compare this sce- 
nario with the smooth background of Figure 1161 but here 
we have to be careful about what quantities we want to 
compare. The recombination rate in a fully ionized region is 
Qsnfj = 1.48 x 10~ 21 s _1 cm~ 3 , so if a fraction /i on of the 
universe is contained in ionized bubbles, then the average 
photon emission coefficient needed to maintain equilibrium 
is j — oib rift /i on . Equating this to the coefficient for the 



© 2007 RAS, MNRAS 000, 000-000 



14 R. H. Kramer and Z. Haiman 



smooth -Xhii, bg = 0.75 case we obtain / ion = 0.75 2 = 0.56. 
In fact we will need an even higher background in the bub- 
bles in order for them to be highly ionized, as we will discuss 
below. 

Next we need to determ ine the size of the ionized bub- 
bles. iFurlanetto et alj l|2004 ) calculated that when 56% of 
the IGM was contained in galactic H II regions, the H II 
region radius distribution peaked at ~ 4 Mpc comoving, or 
about 0.6 physical Mpc at z = 6. We have set up a back- 
ground ionization profile consisting of evenly-spaced ionized 
regions with a diameter of 1 Mpc, occupying 56% of the 
IGM volume (and therefore separated by neutral regions of 
0.79 Mpc along the LOS). The neutral fraction within the 
bubbles is initially set to 10~ 3 . The background needed to 
maintain this level of ionization is 



3 



(To r (1 — Xhii, bg) 



2.01 x 10" 21 s _1 cm" 3 



(19) 



where r = 0.5 Mpc is the radius of the bubble. This is 35% 
higher than the background in the smooth case. We again 
use equation [15] to approximate radiation transfer for the 
background photons. 

The resulting ionization structure is shown in Figure 
1171 Without the background, the fronts in the inter-bubble 
regions would be essentially the same shape as they are in to- 
tally neutral IGM, except sliced and shifted apart where the 
bubbles occur, since the ionizing radiation traverses the bub- 
bles unabsorbed (except by helium and the small amounts 
of residual hydrogen) so the front can pick up where it left 
off on the other side. The presence of helium and the ge- 
ometric dilution cause some changes in the shape, but the 
background radiation transfer has a greater effect. Once the 
front has swept past a bubble, it's background photons can 
escape, and they stream out to contribute to the ionization 
of the inter-bubble regions. 

This type of background structure may have interest- 
ing observational manifestations. Within the inter-bubble 
regions, the front is somewhat thinner than with either no 
background or a smooth background. The front is thinner 
than in the no-background case because the background 
photons contributing to the movement of the front have 
a short MFP. The front is thinner than in the smooth- 
background case because it is "trapped" by the fully neutral 
inter-bubble regions where the mean free path is short, so 
it is slowed down, meaning that a more highly ionized part 
of the front (e.g. r_3) can "keep up" better with a more 
neutral part (e.g. r_2.s), resulting in a thinner front over 
all. There will be brief windows when r_3 is still trapped 
in the residual neutral gas on one side of a bubble, while 
r_2.5 is on the far side, but it would require detailed simula- 
tions of absorption spectra to determine the observability of 
such a situation. The shape of the front at neutral fractions 
more relevant to 21 cm observations again tends to suggest 
smaller thickness measurements, but it seems likely that the 
thicknesses could still be measurable. 

The signature of the hard quasar spectrum is also still 
quite visible on large scales, for instance in the separation 
between Xhi ~ 10 -4 and Xm ~ 0.1. This suggests that a 
combination of 21 cm and Lyman-series observations, par- 
ticularly when smoothed on large scales or averaged over 
many quasars, could place constraints on the ionizing SED 
of the quasar or quasars. If the patchy ionization structure is 
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Figure 18. Comparison of the neutral fraction with a uniform 
ionizing background (solid) and a biased background from pre- 
existing galaxies clustered around the quasar (dashed). Approxi- 
mate background RT is used in both cases. Spectral parameters 
are the fiducial ones used in FigurefS] The three sets of curves cor- 
respond to source ages of log(f/years) = 4, 7 and 8, respectively. 
In the first panel, the central pre-ionized bubble has a radius of 
1 Mpc, while in the second panel the radius is 2 Mpc. 



resolved in future, high-resolution 21cm observations with 
SKA, then the thickness of the front in-between the indi- 
vidual galaxy-bubbles can reveal evidence for a hard quasar 
spectrum separately along each direction measured from the 
quasar. 

Finally, we consider the effect of clustering o f back- 
ground sources around the quasar. IWvithe fc Loebl (I2007T ) 
have explored the mean ionization structure around a quasar 
resulting from the biased clustering of galaxies around its 
host halo (not addressing the discreteness of the bubbles in 
the radial profile). They point out that even outside the I- 
front, the IGM is ionized 10-20% more than the global mean 
ionized fraction. While a high-resolution 21cm interferom- 
eter, such as SKA, could resolve such clusters of galaxy- 
bubbles outside the quasar's H II region, the next genera- 
tion of experiments will likely not have sufficient resolution. 
In this case, the I-front will appear extended due to the 
presence of galaxies. 

In order to explore the effect of such a pre-ionization 
profile on our conclusions, we implemented an ionizing back- 
ground designed to qualitative ly mimic the backgroun d ion- 
ization structure predicted bv lWvithe fc Loebl l|2007h when 
the IGM far from the quasar is at Xhii, bg ~ 0.7-0.8. In 
this scenario pre-existing galaxies produce an H n region 
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centered on the quasar; outside this H n region, there is 
a tail of excess partial ionization, slowly approaching the 
background level of the neutral fraction in the IGM. 

To produce the first panel in Figure 1181 we set up an 
ionizing background to produce a 1 Mpc radius inner ionized 
region (initially set at Xhii = 2 x 1CF 3 with j = 4.41 x 



10 21 s 'cm 3 ). Outside of that, the profile is set to 



Xmi,BG = 0.75 + 



(r/Mpc) + 1.91 



1.81 



(20) 



with j — as n n X H n BG for r > 1 Mpc. In the second panel 
the radius of the pre-ionized region is increased to 2 Mpc 
and the outer profile is shifted accordingly. 

Figure [18] illustrates what happens to the I-front as it 
passes through such a biased background. The dashed curves 
show the I-front expanding into the pre-ionized region, con- 
trasted with the case when the pre-ionization is uniform 
(solid curves). The figure reveals that once the front has 
passed beyond the high-background inner region, it quickly 
begins to approach the shape of the ionization profile in a 
uniform background, though the extra flux from the sources 
clustered close to the quasar keeps it ahead of the front 
in the uniform case. The presence of helium is particularly 
important in this case because it determines the speed at 
which the front can propagate through the pre-ionized re- 
gion (since helium is assumed to be unaffected by the back- 
ground and therefore still neutral). We have actually cho- 
sen rather sma ll sizes for our central pre-ionized region. 
IWvithe fc Loebl (|2007f ) predict central regions of r ~ 2-8 
Mpc when Ami, bg ~ 0.7-0.8. A larger region implies a 
larger background flux once the front has passed the inner 
region. This would prolong the time during which the front 
shape is affected by the background bias. Even a region as 
small as 2 Mpc results in fronts whose inner parts are signif- 
icantly less broad at 10 s years, as seen in the second panel. 



5 CAVEATS 

Before we summarize our conclusions, let us briefly review 
the limitations of the methods we have employed in reaching 
them. 

The one-the-spot approximation — in which we assume 
that ionizing photons produced by recombination are ab- 
sorbed at their point of origin — is strictly valid only where 
the mean free path of the recombination photons is much 
shorter than other length scales of interest, such as the thick- 
ness of the front This condition is violated where the neu- 
tral fraction is very low, in the inner part of the ionized 
region. It is also less accurate for helium- derived photons, 
since they are generally more energetic. The result is a slight 
overestimate of the ionization rates, and therefore a slight 
underestimate of the neutral fraction, close to the quasar. 
The effect is less important for harder spectra (because they 
leave higher neutral fractions overall), and is unimportant 



5 The situation is slightly different for soft (e.g. stellar) spectra, 
in which the mean ionizing photon from the source is also close 
to the threshold. In that case both the source spectrum and the 
recombination radiation are concentrated near the threshold. See 
lOsterbrock fc Ferlanj J2006T) . 



in the outer part of the front. For instance, the mean free 
path at 13.6 eV, Zmfp(^o) < 0.01 Mpc for Xm > 0.07. 

We have ignored the propagation time for photons (due 
to the finite speed of light). Our results accurately predict 
the observed structure along a line of sight aimed directly at 
the quasar, in an absorption spectrum, for example. To sim- 
ulate two dimensional projections of the ionization structure 
(like 21 cm maps), we would need to apply a simple trans- 
formation (see citations in § 12.111 . In fact, since the thickness 
of the I-front varies with time, in principle, a comparison of 
the thickness in transverse and radial directions, in addition 
to the angular dependence of the overall bubble size due 
to finite-speed-of li ght effects (|Wvithe et al.ll2005l ; lYull2005l ; 
IShapiro et alB ooei, will contain information about the age 
of the quasar. 

We have assumed the density of the IGM surround- 
ing the quasar is homogeneous on a macroscopic scale, with 
smaller variations accounted from by a uniform clumping 
factor. Quasars are actually expected to form in biased en- 
vironmental conditions, with high gas density and a high 
density of neighboring halos. The largest effect would be 
along lines of sight that intersect high-density Lyman limit 
systems, which could stall the propagation of the I-front and 
shield the gas further along the LOS from ionizing photons. 
The less extreme effects of density inhomogeneities would be 
somewhat ameliorated by certain observational necessities. 
In order to get a large enough signal to noise ratio to mean- 
ingfully constrain the structure of quasar H II regions, it may 
be necessary to stack 21 cm images of multiple bubbles. Sim- 
ilarly, measurements of the spectra of many quasars would 
be combined in order to constrain the ionization structure 
with Lyman-series absorption observations. In either case, 
density inhomogeneities would add to the random "noise", 
but could, to some extent, be averaged out. Width measure- 
ments may be more robust under suc h conditions than th e 
simple size measurements discussed bv lMaselli et~aH (|2007h . 
however more study is needed to determine that conclu- 
sively. 

As mentioned in section [2~4l we have assumed the IGM 
is isothermal, justifying this by pointing out that most of 
the heating will occur before recombinations become impor- 
tant. Therefore, as long as we choose the gas temperature 
to match conditions expected inside the H II region, this 
assumption should have little effect on our results. 

Finally, we should remind the reader that since we did 
not actually simulate absorption spectra, we have not taken 
into account the effect of the GP trough damping wing (from 
the neutral hydrogen surrounding the ionized region), which 
could be very impor tant for the observational de termination 
of Rp/R a . In fact, Me singer fc Haimanl (|2004T ) essentially 
use this effect to constrain the neutral hydrogen fraction in 
the IGM. 



6 DISCUSSION AND CONCLUSIONS 

We have found that high-redshift quasars residing in a par- 
tially neutral IGM could produce ionization fronts with an 
observable thickness, given a sufficiently hard ionizing spec- 
trum. We discovered that simulating the time-dependent 
evolution of the front is crucial for making accurate pre- 
dictions of the ionization structure, and that the presence 
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of helium in the IGM and secondary ionizations by high- 
energy photoelectrons can both have significant effects on 
the ionization structure. 

With an intrinsic hydrogen column density 
\og(Nu/cm~ 2 ) > 19.2 or a sufficiently hard power-law 
spectrum combined with some obscuration (e.g. spectral in- 
dex s< 1.2 at log(Nu/cm~ 2 ) > 18.0), the outer thickness of 
the front exceeds ~ 1 physical Mpc and may be measurable 
from the three-dimensional morphology of its redshifted 
21cm signal. 

The highly ionized inner part of the front, which 
may be probed by Lyman line absorption spectra, remains 
thin for bright quasars unless a large obscuring column 
(log(iVH/cm -2 ) > 19.2) removes most of their ionizing pho- 
tons up to w 40 eV. Otherwise the lowest energy photons 
with the shortest mean free paths always dominate ioniza- 
tion at the inner face of the front resulting in a small thick- 
ness. 

Highly absorbed ionizing spectra leave a relatively large 
neutral fraction within the H II region, which means that 
the Lyman-series optical depths can be large even within 
the front. For sources with log(iVH/cm _2 ) > 19.8, the Ly- 
man a trough (where the neutral fraction is > 10 -3 ) un- 
derestimates the size of the H II region by a factor of > 4. 
The bias can get as large as a factor of ~ 8 within out pa- 
rameter space. This is in addition to other effects that bias 
this measurement by a muc h smaller amount, as already 
discussed by several authors llMaselli et al]|2007l ; iLidz et all 
l2007l ; lBolton fc Haehneltll2007al ). 

These obscured spectra also result in a large difference 
between the sizes measured by the Lyman a trough, and 
those measured by the Lyman /3 trough (where the neu- 
tral fraction is >10~ 2,5 ). We estimate that Rp/R a could 
be higher than 1.2 for log(iVH/cm _2 ) > 19.8 with a quasar 
age of 3 x 10 7 years. 

We explored the effects of uniform and non-uniform ion- 
izing backgrounds, finding that even with a large uniform 
background the thickness of the ionization front has the po- 
tential to constrain the spectral parameters of the quasar. 
An ionization front propagating into a "swiss-cheese" back- 
ground of small pre-ionized bubbles obviously makes for a 
more complicated ionization structure, but again, the shape 
of the front on both large and small scales is controlled 
largely by the quasar's ionizing SED. The greatest chal- 
lenge is presented by a quasar turning on within a large 
pre-existing ionized region due to galaxies clustered around 
the quasar. In this case, the thickness of the front does not 
reflect the quasar's spectrum until the quasar has pushed 
the front a significant distance farther into the IGM. A mea- 
surement of a thin front around a quasar cannot, therefore, 
be interpreted as definitive evidence of a soft or unabsorbed 
quasar spectrum. On the other hand, measuring a thick front 
does suggest a highly absorbed quasar spectrum, because the 
galactic sources producing the pre-quasar ionization should 
be soft. 

IShapiro et ail l|2004h and iThomas fc Zaroubil (j2007T l 
have also used 1-D radiative transfer calculations (with 
slightly different numerical methods) to explore high- 
redshift ionization structures, with t he addition of self- 
consis tent temperature calculations. IThomas fc Zaroubil 
also included collisional ionization, and an evolving 



mean IGM density, while IShapiro et ail ((2004) (who wore 
focused on the evaporation of minihalos by the I-front) in- 
clude d gas dynamics and trac e metals. Like the present 
work, |Th omas fc Zaroubil l|2007t) included secondary ioniza- 
tions, and both papers included helium, and calculations of 
the time evolution of the front. Neither paper included the 
ionization of hydrogen by helium recombination radiation, 
and they explored only blackbody (stellar) and power-law 
or truncated power-law (quasar or "miniquasar") sources, 
rather that the ultra luminous obscured quasars that we are 
interested in here (though the miniquasar spectra ar e fairly 
similar to some of our spectra). IShapiro et al.l l|2004l ) found, 
as we did here, that I-fronts around sources with harder 
spectra are thicker. However, they studied sources at higher 
redshift (z — 9), which were much less luminous, and the 
hardest spectrum they examined (a 10 5 K blackbody, rep- 
resenting pop III stars) was still softer than most of the 
spectra in the range included in the present paper. As a re- 
sult, they found fronts much thinner (up to ~ 0.01 Mpc; 
see their Figures 7 and 8) than most of those discussed in 
the present paper. Amon g the many interesting findings of 
IThomas fc Zaroubil (|2007f ) are several relevant to the present 
study. Even with their harder power-law spectrum (s = 1 
extending from 200-10 4 eV), they did not find extended ion- 
ization tails outside of the front, due to the combination of 
shorter lifetimes and higher redshift (when the IGM density 
was higher). They do, however, find an extended kinetic tem- 
perature structure coupled to the spin temperature so as to 
be observable with 21 cm instruments, and find that such 
observations could be useful in discriminating between var- 
ious source spectra. 

Finally, we found that the contours of degeneracy in our 
parameter space for measurements using 21 cm observations 
are oriented differently from those for measurements using 
Lyman-series absorption spectra. This suggests that studies 
combining both types of observations have the potential to 
break the degeneracy and constrain both the intrinsic ab- 
sorption column and spectral index of a quasar's ionizing 
radiation at z > 6. 
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